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Abstract 

Three-loop quantum corrections to the effective action are calculated for 
' = 1 supersymmetric electrodynamics, regularized by higher derivatives. 
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Using the obtained results we investigate the anomaly puzzle in the considered 
model. 



1 Introduction. 

^ ! It is well known [1, 2, 3, 4], that in supersymmetric theories the axial and the trace 

^ I of the energy-momentum tensor anomalies are components of a chiral scalar super- 

multiplet. Adler-Bardeen theorem [5, 6] asserts that there are no radiative corrections 
to the axial anomaly beyond the one-loop approximation, while the trace anomaly is 
proportional to the /5-function [7] to all orders. Therefore it seems to imply, that the 
/3-function in supersymmetric theories should be exhausted by the first loop [8]. It does 
take place in models with N = 2 supersymmetry [9]. However explicit perturbative 
calculations find higher order corrections to the /3-functions of = 1 supersymmet- 
ric theories, regularized by dimensional reduction [10, 11, 12]. This contradiction is 
usually called " the anomaly puzzle" . 

Many papers were written in the attempt of solving the anomaly puzzle in su- 
persymmetric theories. For example, in [13] the anomaly puzzle is argued to be a 
consequence of the difference between the usual and Wilsonian effective actions. In 
particular, the authors noted, that there was a nontrivial contribution to the /3-function 
related with the Konishi anomaly [14, 15]. The investigation of this contribution in 
[13] and the investigation of instanton contributions in [16] have led to the construction 
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of the so-called exact Novikov, Sliifman, Vainshtein and Zakharov (NSVZ) /3-function. 
For N = 1 supersymmetric electrodynamics (SUSY QED) considered in this paper the 
NSVZ /3-function has the following form: 

/3(a) = ^(l-7(«)), (1) 

where 7(0;) is the anomalous dimension of the matter superfield. Explicit perturbative 
calculations with the dimensional reduction (DRED) verify the NVSZ /^-function up to 
the two-loop order. Nevertheless, the three-loop results obtained in [17, 18, 19] do not 
agree with the NSVZ /3-function. However [18] this disagreement can be ehminated by 
a special choice of renormalization scheme, the possibility of such a choice being highly 
nontrivial [20]. Actually it is possible to relate DRED scheme and NSVZ scheme order 
by order [21] in the perturbation theory. It is worth mentioning, that at two-loops the 
NSVZ /3-function was also obtained with differential renormahzation [22] . For example, 
for = 1 SUSY Yang-Mills the calculation was made in [23] . 

However the relation between F and the Wilsonian action remained unclear. This 
problem was avoided in another solution of the anomaly puzzle, proposed in [24] . The 
main idea of [24] is that the higher order corrections in NSVZ /9-function are due to 
anomalous Jacobian under the rescaling of the fields done in passing from holomorphic 
to canonical normalization. In the case of supersymmetric electrodynamics holomor- 
phic normalization means, that the renormalized action is written as 

Sren = ^^3(A//i)Re J d^X d'^O Wa^'^Wf, + 

+Z{A/i^) ^ J d^x d^9 (0*6^^0 + 0*e-2^0) , (2) 
while in the canonical normalization 

Sren = ^^3(A///)Re j d'^X d^O WaC'Wh + 

+ ljd''x d^e {(l)*e^^4> + 4>*e-^^4>) . (3) 

In the former case the /3-function is supposed to be exhausted at the one loop, while in 
the latter one it coincides with the NVSZ result. In principle this solution is different 
from the one, given by Shifman and Vainshtein. Moreover, it contradicts the results of 
exphcit two- loop calculations, made with DRED. 

It would be natural to suppose, that in the holomorphic normalization the /3- 
function is exhausted at the one-loop if higher covariant derivative regularization 
[25, 26], supplemented by the Pauli-Villars, is used. This regularization is known to 
yield the same result for one-loop logarithmic divergences as the dimensional regular- 
ization (or dimensional reduction) [27] . The explicit two-loop calculations for theories, 
regularized by higher derivatives (HD), were made first in [28, 29] for A" = 1 SUSY 
QED and ga\x> a zero tiA o-loo]) contribution to the /3-function defined by ^ 

^Note, that it is the /3-function, that is proportional to the trace anomaly. 
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This result implies the absence of the anomaly puzzle in view of the solution pro- 
posed in [24] . However it was not quite clear why different regularizations give different 
results for the scheme independent two-loop /3-function. Actually in [28] wc noted, that 
the using of the HD regularization leads to a nontrivial contribution of diagrams with 
insertions of one-loop counterterms, which does not exist for the dimensional reduc- 
tion. The calculations of this contribution with different regularizations were analysed 
in [30], where the difference of the results for the scheme independent two-loop (3- 
function was attributed to the mathematical inconsistency of DRED [31], which had 
been pointed in [32]. In particular, the inconsistency of DRED leads to incorrect zero 
results for anomalies, because DRED does not break the chiral symmetry. It is neces- 
sary to stress an essential difference between dimensional regularization (DREG) [33] 
and DRED: DREG allows to derive the axial anomaly unambiguously [33]. However 
DREG explicitly breaks supersymmetry and is not convenient for the calculations in 
supersymmetric theories. Let us note, that anomalies can in principle be calculated 
with DRED. However for this purpose it is necessary to impose mathematically incon- 
sistent conditions like ti:{AB) ^ ii{BA) [34] or use some identities between 7-matrices, 
which are vahd only for n > 4 [35]. (DRED requires that the space-time dimension 
n should be less than 4 [31].) However such conditions can not be imposed if the 
calculations are made by the supergraph technique. Hence the axial anomaly and the 
Konishi anomaly, calculated with DRED, are equal to 0. As a consequence the addi- 
tional anomalous contribution, pointed in [13], is omitted if the theory is regularized by 
DRED. HD regularization is mathematically consistent and allows to calculate anoma- 
lies correctly. In particular, the anomalous contribution to the /3-function, obtained 
with HD, is not equal to 0. Actually this contribution is a sum of Feynman diagrams 
with insersions of counterterms on matter lines. The sum of such diagrams is equal to 
with DRED and agrees with the results of [13] and [24] with HD regularization. After 
rescaling, which converts (2) into (3), the diagrams with insersions of counterterms 
vanish, and the /^-function becomes equal to the NSVZ expression. 

It is necessary to note, that although the /9-function (4) is exausted by the first 
loop in the holomorphic normalization, the Gell-Mann-Low function has contributions 
from all orders. This contradiction is discussed in the present paper. We argue, that 
if the Adler-Bardeen theorem is valid and the bare coupling constant does not depend 
on /i, then the generating functional depends on ^ due to the rescaling anomaly and 
/9-function (4) is not related with the Gell-Mann-Low function. Therefore there is 
no contradiction between the form of the Gell-Mann-Low function and the multiplet 
structure of anomalies. 

One more purpose of this paper is the calculation of the /3-function in the three- 
loop approximation. It is desirable in order to avoid some possible errors or incorrect 
interpretation of the results, especially if we take into account, that the three-loop 
/9-function, considered as a function of a, is scheme-dependent. The thrcc-Ioop con- 
tribution to /5-function (4) is found to be 0, and agrees with the predictions of [24] 
and [30] . It is worth mentioning, that in the three-loop approximation the sum of the 
diagrams without insersions of counterterms (on matter lines) for a large number of 
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subtraction schemes is equal to the exact /9-function (calculation with DRED gives the 
NSVZ /3-function only after a redefinition of the coupling constant). The sum of the 
diagrams with insersions of counterterms in two- and three-loop approximations agrees 
with the exact expression found in [30], and cancels the other two- and three- loop 
contributions. 

The paper is organized as follows: 

In section 2 we consider = 1 SUSY QED and regularize it by higher deriva- 
tives. The three-loop /3-function and its relation with two-loop anomalous dimension 
are analysed in section 3. In particular, the three-loop contribution to the /3-function 
is found to be 0. In section 4 we explain why the results are different from those ob- 
tained with DRED. The anomaly puzzle is considered in section 5. Section 6 contains 
some concluding remarks. The details of the calculations are presented in appendixes. 
Appendix A contains expressions for various groups of Fcynman diagrams. The calcu- 
lations of the corresponding contributions are made in appendix B and the most useful 
three-loop integrals are analysed in appendix C. 

2 N = 1 supersymmetric electrodynamics and 
higher derivative regularization. 

= 1 supersymmetric electrodynamics is described by the following action: 

So = ^Re J d*x d^e WaC^Wb + \j d^x d^9 (0*6^^0 + cf>*e-^^^) . (5) 
Here and (p are chiral superfields 

4>iy, 9) = ifiy) + ^(1 + 75)^(2/) + lo{l + 75)^/(2/); 

9) = <p{y) + 9{1 + ^^iv) + ^^"(1 + 1^)9 f{y), (6) 
where y^ — + i9^'^^^9 /2. Two Majorana spinors -0 and -0 form one Dirac spinor 

* = ^((l + 75)V^ + (l-75)^)- (7) 
V in (5) is a real superfield 

V{x,9) = C{x)+iV29^,i{x) + \{99)K{x) + '-{h^9)H{x) + ]^(9Yl^e)A^{x) + 

+V2(99)9{n,x(x) + Iri^d.^x)) + \{09r{D{x) - \d'C{x)), (8) 

where, in particular, is an Abelian gauge field. The superfield Wa in the Abelian 
case is defined by 
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Wa = -D{1 - 75)L'[(1 + ^^)DaV 



(9) 



where 



(10) 



is a supersymmetric covariant derivative. 

In order to regularize model (5) by HD its action should be modified as follows: 



5*0 — >■ 5* — 5*0 + S'a — 

= —ReJd^xd^0W.C^'{l + —)w,+ 

+ \j d^x d^9 (0*e2^0 + 0*e-2^0) . (11) 

Note, that in the Abehan case the superfield is gauge invariant, so the higher 
derivative term contains usual derivatives. 

The quantization of (11) can be made by using standard technique described in [36] 
and is not considered here. It only needs mentioning that the gauge invariance was 
fixed by adding 

where 

D' = ^D{1 + ^,)D; D' = ^D{1-^,)D. (13) 

After adding such terms the free part of the action for the superfield V is written in 
the simplest form 

S„ + Sgf^^l d'x d'e Vd' (l + ^) V. (14) 

In the Abelian case diagrams containing ghost loops are missing. 

The superficial degree of divergence for the model (11) is (see e.f. [28]) 

a;A = 2-2n(L-l)-E^(n + l), (15) 

where L is a number of loops and is a number of external 0- lines. According to (15) 
divergences remain in one-loop diagrams even for n > 2. In order to regularize these 
divergences it is necessary to insert Pauli-Villars determinants [6] into the generating 
functional. Due to the supersymmetric gauge invariance 

V^V-^-{A + A-^)- cP^eU; 4>^e-% (16) 
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where A is an arbitrary chiral scalar superfield, the renormahzed action can be written 
as 



Sren = ^^3(A/M)Re J X (f 6 W aC''' { 



Q2n, 
A 2r) ' " 



+Z{K/^i) ^ j d^x d^e {(f>*e''^(f> + 0*e- 



-2V 



(17) 



Hence the generating functional is 



f ~ ell/" (7 

j DV D(l)D^l[(^detPV{V,Mi))\x.pY j-^ j d'^xd'^OVd'^il + j^)v - 
-4i^(^3(A^) - l) / d^xd^9VIiy,d\l + —)V + 



+ ^Z(A//i) I d^xd^eU*e^^(l) + 4>*e-^^ 



+ 



+ 



J d^x dH JV + J d^x d'^e (j (f) + j$) + J d^x d^e {j*(f)* + ~f4)*) 



(18) 



where 



(detPF(y,M)) ^ = jD^D^expli Z(A//i)i ^ rf^xrf^^ ($*e^^$ + 



+ d^xdHM^^ + d'^xd^eM^*^* 
and the coefficients q satisfy equations 



(19) 



(20) 



Below we assume, that Mi — a^A, where Oj are some constants. The insertion of 
Pauli-Villars determinants allows us to cancel remaining divergences in all one-loop 
diagrams, including diagrams with insertions of countcrtcrms. Later we will show, 
that the divergencies in the sum of two- and three-loop diagrams with Pauli-Villars 
loops cancel each other. Therefore, for diagrams with loops of Pauli-Villars fields it is 
unnecessary to introduce any other regularization. 

In our notations the generating functional for connected Green functions is defined 

by 



W = -i\nZ, 

and an effective action is obtained by making a Legendre transformation: 



(21) 



V - j d^xd'^eJV - j d^x d^9 (j (j> + j4>)- j d'^x d^e + , (22) 
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where J, j and j is to be ehminated in terms of V, (f) and 0, through solving equations 

5W_ 5W ~ 5W 

^=1/' "^^W ^ ^ 

After obtaining Sren, it is possible to find the /3- function and the anomalous dimen- 
sion, which in our notations are defined by 

^ _ d /en _ d\nZ 

(We assume, that the bare coupling constant Cq, defined by 

i = 4^3(A^), (25) 

does not depend on /j,. Hence the rcnormalizcd coupling constant e depends on /x.) It 
is easy to see [7], that the trace anomaly is proportional to /3-function (24). 

The /3- function and anomalous dimension, given by (24), which are considered as 
functions of e, are changed at the simultanious redefinition of the renormalized coupling 
e and the renormalization constant Z^, provided = const. In other words they 
depend on the renormalization scheme. If f3 and 7 are expanded in powers of e^, then 
the coefficients of the /3-function and anomalous dimension become scheme-dependent 
starting from the three- and two-loop approximation respectively. 

Note, that it is possible to use another definition of the /3-function. Let us consider 
transversal part of the two-point Green function for the gauge field: 



XT/ XT/ (^Pm^" + ^Qi^y^) = 

^{27ry6\p + q)p'U,/,5\e,-ey)d-\a,i,/p), (26) 



where 

1 



ni/2 ^ -j^,D^D'Ca,DK (27) 
Then it is possible to define Gell-Mann-Low function 

d 

p(d{a,x)j = —x—d{a,x). (28) 

Taking into account, that the effective action should not depend on the normalization 
point /I and differentiating equation (26) over In /i we obtain 

d 

= (3{d{ a,x)j—(3{a)—d{a,x). (29) 



In particular at x = 1 we have 



(3{a) ^ (3{a)—, (30) 
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where a = d{a, 1). Therefore, if the generating functional does not depend on /i, then 
both definitions of the /^-function are equivalent. 

In order to find the /3-function and the anomalous dimension it is necessary to 
calculate all IPI graphs in the considered approximation. The expressions for them 
are constructed in accordance with Feynman rules, which can be formulated as follows: 

1. External lines give the integration 

n / ''(^.v) / '^*K) Ep.). (31) 



E 



where E runs over external momenta. 
2. Propagator of the superfield V is 



^""^ S'iOi-e^). (32) 



(A;2 + i0)(l + (-1)"A;2VA2") 

3. Massless (j) — (f)* and ^ — 0* propagators are 

77^ rD^D^S^Oi - 62) . (33) 

16(A;2 + iO) V 1 y ^ 

(Note that the considered action is quadratic in matter superfields and Feynman rules 
can be simplified in comparison with, say, Wess-Zumino model.) 

4. Pauli-Villars fields are present only in the closed loops. Each internal line $ — 
or l> — l>* corresponds to 

^ D''DH\di-d2), (34) 



16(A;2 - M2 + iO) 
and each internal fine $ — l» or $* — l»* - to 

4(P-1^ + «) ^'^'W-^-) and ,(,._^^,o) g-^'W-fe) (35) 

respectively. For each loop of Pauli-Villars fields it is necessary to introduce — q. 



{2ny 



5. Each loop yields an integration over a loop momentum j 

6. Each vertex gives J d'^O. 

7. There are usual combinatoric factors, which can be found from the generating 
functional (18). 



3 Three-loop /^-function. 



The diagrams contributing to the three-loop /3-function are presented in Fig. 1 - 
8. Note, that each internal matter loop in these diagrams can correspond to and 
(p fields or to Pauli-Villars fields. We devided Feynman diagrams into some groups 
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and presented the expressions for all these groups in appendix A. Then the three-loop 
correction to the effective action can be written as 



ArP = Re / d'e ^^WMC^'W.i-p) X 



(27r)^ 

X i^fi-ioop + /o + fopv + foo + foopv + foo + fxx + fx + fx2j , (36) 

where 

fi-ioop is a contribution of one-loop diagrams, presented in Fig. 1; 

fo is a contribution of diagrams, presented in Fig. 5, containing a single loop of 
the superfields and 0; 

fopv is a contribution of the same diagrams, having a loop of Pauli-Villars fields; 

foo is a contribution of the two-loop diagrams, presented in Fig. 2 (without Pauli- 
Villars fields) and the three-loop diagrams, presented in Fig. 4, which contain an 
internal loop of matter superfields or an insersion of one-loop counterterms on the 
photon line and external loop of and superfields; 

foopv is ^ contribution of the same diagrams with the external loop of Pauli-Villars 
fields. 

foo is a contribution of the other diagrams with two loops of the matter superfields, 
presented in Fig. 3; 

fx is a contribution of diagrams with an insersion of two-loop counterterms, pre- 
sented in Fig. 6; 

fxx is a contribution of diagrams with two insersions of one-loop counterterms, 
presented in Fig. 8; 

fx2 is a contribution of diagrams with an insersion of one-loop counterterms, pre- 
sented in Fig. 7. 

The explicit expressions for all these contributions were found by means of calculat- 
ing of the corresponding Feynman diagrams. To check correctness of these calculations 
we verified the cancellation of noninvariant terms, proportional to 

/ -0j-,d"(^ViP^(^)n-P,O). (37) 

The results, presented in appendix A, are analysed in appendix B. Let us briefiy discuss 
them: 

1- foo = 0, because the substitution in a loop changes the sign of a diagram. 
Indeed, in this case the diagrams, having the same superfields in both loops (0 and (p 
or (f) and 0), are cancelled by diagrams with loops of different superfields {(f) and 0). 
For the diagrams with Pauli-Villars fields the result is the same, but its derivation is 

more complicated. 

2. The sum of fxx, fx and fx2 agrees with the exact expression for the sum of 
diagrams with insersions of counterterms 

- In Z Re / d'^x d^O WaC^Wb + finite terms, (38) 

found in [30]. The corresponding contribution to /3-function in the considered approx- 
imation is 
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A/3 =—7(0;). (39) 

TT 

According to the results of the one-loop calculations and the predictions of the 
renormgroup (see e.f. [28]) the constant Z is given by 

Z(A//i) = l + -(ln-+5i) + 

+ ^ ( In^ - + ^1 In -) - 72 «^ ( In - + + 0{o?). (40) 

Here 720;^ is a two-loop contribution to the anomalous dimension and we assume, that 
at the one-loop the counterterms are 

+^(ln^ + ^i) / d^xd^^(0V^0 + 0*e-^^0), (41) 



where bi, gi and g2 are arbitrary finite constants, which define a subtraction scheme. ^ 
3- foopv s-nd fopv finite and do not contribute to the divergent part of the 
effective action. This means, that the sum of all diagrams with Pauli-Villars loops is 
finite, although there are divergences in some of such graphs. However, Pauli-Villars 
regularization always assumes the existance of divergent diagrams and the cancellation 
of the divergences between different graphs. Therefore, in the considered case it is not 
necessary to introduce any more regularization. 

4. The analysis of foo and fo is rather involved, because the corresponding integrals 
are very complicated. Each of these integrals depends on A/p and is the sum of a third 
degree polynomial in In A/p and a function, finite at A — > 00 (or equivalently at p — > 0) 

/(A/p) = /3 In^ - + /2 In^ - + /i In - + /o + o(l). (42) 
p p p 

Let us assume, that the limit 

(43) 



A^/(A/rt 



p=0 

exists. Then /s = and /2 = 0, while the considered limit is /i. 

In appendix B we prove, that for foo and fo the limit (43) exists and 



Oo 



o 



p=0 



167r2 



(44) 



p=0 



^It is convenient to include the higher derivative term in counterterms, because its presence sim- 
pUfies an expression for the two-loop anomalous dimension [29] . In principle, this term is not essential 
and can be omitted. 
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where 



^2 



+ 



J (2^0" (2^ + pf{k + Z + A;27A2") (l + f "/A^") 



+ 



r d K d I 

J (27r)4 (27r)4 ^^2^^. + p)2/2(/ + p)2(^ + / + py(l + k^n//^2n^ (l + /2n/A2n^ 

- / ^'''^ J [ ^ - (46) 

- E Q / 7^ r^^ V - ^ ( In - + fei) (l + ^t^VA^")! . 

^ 'J {271^ ^12 + Mi)({k + ir + Mi) Stt^V / J) 

Prom (44), (45) and (46) we see, that the integral over three loop momenta is reduced 
to the integral over two loop momenta. It is very nontrivial, that can be seen from 
the calculations, done in appendixes B and C. In our opinion these facts confirm the 
correctness of the obtained results. 

Note, that Zi and Z2 are present in the two-loop two-point Green function for the 
matter superfield [29]: 



(2ny 

2^1 ( In h Oi ) -|- (In h Oi ) -|- two-loop counterterms > . (47) 



47r2 

This expression can be formally written as 

r^'^ = i / (04 {f^v. e) H-p, e) + 0*(p, e) ^{-p, e)) x 

xl^l- zi- Z2- Azizi + Azi + Az2 + A(Azi ^i) |, (48) 

where the operator A is constructed as follows: If / is a function oi A/p and A/fi, then 
by definition A/ is a counterterm, which cancels a divergence of the function /. For 
example, 

Azi^-(\n- + gA. (49) 
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Below we will assume, that the operator A is linear. (In a general case this operator 
can be nonlinear). 

From (48) we see, that the two-loop renormalization constant for the matter super- 
field is given by 

Z ^l + Azi + Az2 + A(^ziZi) + 0(a'^), (50) 

so that 

InZ = Azi + Az2 + A{Azi zi) - ^{Azif + 0{a^) = Azi + Az2 - 
-Ia{z, - Az,y + lA{zf) + lA{{Az,r) - ^-{Azxf + 0{a% (51) 

Taking into account, that due to the definition of A the expressions {zi — A^^i)^ and 
A^(A^i)^j — (A^i)^ are finite, InZ can be presented as 

In Z = A(^Zi + Z2 + 2^i) + finite terms. (52) 

Then the sum of diagrams, defining two- and three-loop contributions to the /^-function 
for A'^ = 1 SUSY QED, for subtraction schemes, corresponding to any linear operator 
A can be written in the following form: 

^"^^^ ^ 16^^^ / d'e ^,W.{p)C^'W,i-p) X 

X + Z2 + -zf — A(^zi + Z2 + 2^i) finite terms^ . (53) 

This expression is finite due to the definition of A, so it is not necessary to add any coun- 
terterms in two- and three-loop approximations. Thus for all renomalization schemes 
with linear A we have: 



^""^ -\n--h + 0{a'). (54) 



Co a(A//i) 



This means, that the two- and three-loop contributions to the /3-function are equal to 
zero and 

p = — + 0{a^). (55) 

71 

Hence the /9-function is exhausted at the one-loop and agrees with the multiplet struc- 
ture of anomalies. 

Note, that the sum of diagrams which do not contain insersions of counterterms on 
matter lines in the considered approximation gives the following contribution to the 
/3-f unction: 
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A/3 = -(l-7(«)). (56) 

This contribution is equal to the NSVZ /3-f unction, but the anomalous dimension is 
cancelled after adding (39), and the final result is comletely defined by the one-loop. 

4 Comparison between HD regularization and 
DRED. 



The /9-function obtained in the previous section is different from the corresponding 
result, found with DRED. In the two- loop approximation the calculations of the efi:ec- 
tive action with DRED and HD were compared in [30]. The difference of the results 
for the /3-function is shown to have originated from the different results for the sum 
of diagrams with insersions of counterterms. With DRED this contribution is 0, while 
with HD it is given by (38). The calculations made in this paper show, that in the 
three-loop approximation we have a similar situation. 

The difference of the results for the sum of diagrams with insersions of counterterms 
[30] is caused by the mathematical inconsistancy of DRED, pointed in [32], because 
this inconsistency leads to zero results for all anomalies. (We assume, that there are no 
assumptions like tr{AB) ^ tr{BA) and all identities are valid for n < 4.) In particular 
the sum of diagrams with insersions of counterterms on the matter lines calculated 
with DRED is 0. 

Let us discuss this in detail: 

In supcrsymmetric theories the axial anomaly is related with the Konishi anomaly 
[14, 15]. Indeed, let us consider 

(57) 

Using equations (6) and (8), it is easy to see, that in components this expression will 
contain (among other terms) 

-^^a^(^7'^75*), (58) 

where the Dirac spinor is defined by (7). It is well known [37], that the conservation 
of the axial current is broken by quantum corrections and in particular 

99 {d^^^r^^"^)) = -09 ^F,,F^\ (59) 

Hence due to the supersymmetry ^ 

''Note, that our arguments can not be considered as a derivation of the Konishi anomaly, because 
(58) does not contain all terms of (57), proportional to 99. A strict derivation of the Konishi anomaly 
can be found in [14, 15]. Our goal is only to remind of the relation between the axial anomaly and 
the Konishi anomaly. 



Im 
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Im {D\<l>*^''(i> + 0*e-2^0)) = ^^ImiWaC'^'W,). (60) 

By performing supersymmetry transformations it is easy to see, that if an imaginary 
part of a chiral superfield is equal to 0, then this superfield is a real constant. Therefore, 
from (60) we obtain, that 



(0*e2^0 + 0*e-2^0) ) = ^WaC^'Wb + const. (61) 



- I d^xD^= I d^xd^e (62) 



Applying 



2 

to (61) and taking a real part of the result, we obtain, that 



Jd^xd^e (0*6^^0 + 4>*e-^^4>)) = -j^Re J d^xd^eWaC'^Wb. (63) 



Because DRED requires, that the space-time dimension n should be less than 4 
[31], it is possible to choose 75 anticommuting with all 7-matrices. Then the chiral 
symmetry is not broken in the regularized theory due to the mathematical inconsistency 
of DRED. As a consequence axial anomaly appears to be 0, while the supersymmetry 
is not broken. Therefore instead of (63) we obtain 



d^x d^e (0*6^^0 + 4>*e-^^4>) ) = 0. (64) 
In DREG such problem can be solved by using 75 with the following properties: 

{75,7m} = 0, ;U = 0,...,3; [75,7^1 = 0, /i > 3. (65) 

Then the chiral symmetry is broken in the regularized theory, and axial anomaly is 
calculated correctly [33]. Nevertheless, DREG breaks the supersymmetry and is not 
well-suited for supersymmetric theories. 

As a consequence of (63) we obtain [30] the identity 



( exp (^i{Z - 1) ^ / d^x d'^e {fe^^4> + 0*6-^^0) j ) = 

= exp ^ - i In Z^^Re j d^x d^O WaC^-^Wb + finite termsj , (66) 

whose l.h.s. is a sum of all diagrams with insersions of counterterms on lines of the 
matter superfield. The corresponding result obtained with DRED, which follows from 
(64), is written as 




Then the sum of all diagrams with insersions of counterterms is 0, that contradicts the 
result for Konishi anomaly. Thus the mathematical inconsistency of DRED gives the 
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result for the /9-function which differs from the corresponding result obtained with HD 
and leads to the anomaly puzzle. 



5 Solution of the anomaly puzzle 



In order to investigate the anomaly puzzle it is convenient to use the higher deriva- 
tive rcgularization, because it is applicable for the calculation of anomalies. 

From the calculation described in section 3 we see, that the /3-function, defined by 
(24) is exausted at the one-loop, while the /5-function defined by (28) has corrections 
from all orders. This seems to contradict equation (30). However, actually there 
is no contradiction, because the generating functional (18) depends on Really, 
due to the rcscaling anomaly (66) it is impossible to remove /x-dependencc by the 
transformation because the anomalous contribution contains /i-dependent 

InZ. Therefore, the /5-functions (24) and (28) are different. (See the derivation of 
equation (30).) The first function is proportional to the trace anomaly and due to the 
multiplet structure of anomalies is exausted by the first loop, while the second one has 
corrections from all orders. 

Nevertheless, it is desirable to avoid /x-dependcnce of the generating functional. It 
can be made by two different ways. We can assume, that the bare coupling constant Cq 
depends on /x. In this case /3-function (24) will have corrections from all orders, but it 
will not be proportional to the trace anomaly. Another possibility is to use canonical 
normalization for the matter superfields and to define the generating functional by 



r ~ 11/" u 

Z = J DV D(j)D4)'[l(^ det 'PV{V, Mi)) ' exp ^ + J^) ^ ~ 

-^(^3(A//.)-l) / d^xd^eVYi,,,d\l + —)v + 



+ j d'^xd^ejv + j d^xd^e(j(i) + ]$) + j d*xd^e(j*(i)* + j*4)*) 

where 



(68) 



( det 'PV{V, M)) ^ = j DI> exp I i d^x d^0 (^^e^^^ + 



d^xd^eMm + =- I d^xd'OM^*^ 



If 



;4™ j2? 



(69) 



Then the /9-function (24) will be proportional to the trace anomaly, but Adler-Bardeen 
theorem is not valid in this case. 
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6 Conclusion 



In this paper we calculated the three-loop /3-function for = 1 SUSY QED regular- 
ized by higher derivatives. Using the standard definition of the generating functional we 
found, that two- and three-loop contributions to the /5-function (24) were for a large 
number of subtraction schemes. In this case the sum of diagrams without inscrsions of 
counterterms on matter lines is exactly equal to the terms of the corresponding order 
in the expansion of the NSVZ /3-function. However two- and three-loop contributions 
are exactly cancelled by diagrams containing inscrsions of counterterms. 

The result for /5-function (24) obtained with HD regularization differs from the 
corresponding result obtained with DRED, because DRED is not mathematically con- 
sistent and does not permit to calculate anomalies [30] (if there are no additional 
assumtions like tT:(AB) ^ ixiBA) e.t.c). In particular, the Konishi anomaly, which 
contributes to /3- function (24), calculated with DRED is 0. This in turn leads to the 
anomaly puzzle. HD regularization enables us to find an anomalous contribution of 
diagrams with inscrsions of the counterterms, which was calculated in [30] exactly to 
all orders and is equal to 

-7(a). (70) 

TT 

The calculations done in this paper confirm this result in the three-loop approximation. 

The result for /3-function (24), obtained with the generating functional (18), is 
consistent with a multiplet structure of the anomalies: Since in supersymmetric theories 
the axial and the trace of the energy-momentum tensor anomalies arc members of a 
supersymmetric multiplet, the /5-function (24) should be exhausted by the first loop. 
In particular, if the theory is regularized by HD the Adler-Bardeen theorem does 
not conflict with supersymmetry, while for theories, regularized by DRED, such a 
contradiction seems to take place [38] . 

However, the generating functional (18) depends on /i due to the rcscaling anomaly. 
As a consequence, the /5-function (28) is different from the one defined by equation 
(24). If we would like to define a /^-independent generating functional, then either 
Adler-Bardeen theorem is not valid or the trace anomaly is not proportional to the /3- 
function. Therefore, if the generating functional does not depend on the normalization 
point, then the arguments based on the multiplet structure of anomalies can not be 
used. In our opinion this solves the anomaly puzzle in the considered model. 

One of the possible ways to define a /i-independent generating functional is the 
using of the canonical normalization (3). Then there are no diagrams with inscrsions 
of counterterms and the /3-function is equal to the NSVZ expression. It is important to 
note, that unlike DRED the HD regularization does not require to tune the subtraction 
scheme. The NSVZ expression (at least in the three-loop approximation) for /9-function 
(24) is authomatically obtained with HD regularization if an operator, constructing a 
counterterm for a given function, is linear. 

It is necessary to note, that so far we considered only the Abelian case. For the su- 
persymmetric Yang-Mills theory the using of higher covariant derivative regularization 
[39] leads to very involved calculations, because in this case Feynman rules become 
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much more complicated. In this case the using of usual derivatives can simplify the 
calculations considerably. However, such regularization breaks the gauge invariance. 
Nevertheless, even in the case of noninvariant regularization it is possible to obtain the 
gauge invariant renormalized effective action by a special choice of subtraction scheme 
[40, 41] . For Abehan supersymmetric theories such scheme was proposed in [42] . Con- 
struction of the invariant renormalization procedure for supersymmetric non-Abelian 
models is in progress. 
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Appendix. 

A Results for Feynman diagrams. 



Having calculated Feynman diagrams, presented in Fig. 1 - 8, we obtained the 
following expressions in the Minkowski space: ^ 
1. One-loop diagrams, presented in Fig.l: 

i r d^k f 1 1 \ 

h-ioo, = - 2 y [k\k+pY ~ ? "\k--Mf){{k+pf-Mf)) ■ ^^^^ 



2. Two-loop diagrams, presented in Fig.2, and three-loop diagrams, presented in 
Fig. 4, with the external loop of and 0: 



X (l + (-l)«A;27A2«)||l + ^(ln^ + 6i)j +2ze2|| j 



dH 1 



dH 



(27r)4 (t2-Mf)({k + tf-Mf), 
3. The same diagrams, with the extenal loop of Pauli-Villars fields: 



(27r)4 t^k + ty 

(72) 



^For simplicity we omit +iO in propagators. 
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X 



(l + (-1)«A;2«/A2«) 
{k + p + qy + — k"^ —p^ 



2 X 



({k + qf - Mf) ({k + p + qf - Mj) (f - Mj) ((g + p)2 - Mj 

+ ^ 



+ 



X 



((A; + g)2 - M2) (g2 _ Mj) ((g + p)^ - Mj) 

f dH 1 



X 



(27r)4 (t2_M2)((A; + t)2-M2 
4. Diagrams with two loops of matter superfields, presented in Fig. 3: 



X' 



■ 2fc2 [{2q + k + lf + {2q + k + l+pf- p") 4^2 



{q + kY{q + lY{q + k+ pY {q + kY{q + k+ p) 
K<l + Pf ^ 4(g + p)^ ^ 4 



2 



(g + A;)2(5 + A;+p)2 + A; +p)2(g + A; + Z +p)2 {k + qY 
2p^ Ak"^ 



{q + kY{q + k + pYiq + k + l+ pY {q + kYil + k + pY 
2kY 2{k + If 



{q + kY{q + k + pY{q + k + I + pY {q + k + I + pY{q + kY 
{k + lfp"^ A{2q + k + l+pY 



+ 



{q + kY{q + k+pY{q + k + l+pY {q + k+pY{q + I) 



+ 



2{k + lY{2q + k + l+pY ^ 2p\2q + k + l + pY 



(73) 



foo = 0. (74) 
5. Diagrams with a single loop of and 0, presented in Fig. 5 



__. r d^k dH d\ 
■'^ ~ ~'J (2^ (2^ (2^ ^ 

k^(l + (-1)"A;2VA2») ^2(^1 + {-l)npn/A.2n^ q2(^q + p)2(^ + A) + 0^ 



(75) 



(q + k + l + pYiq + lY{q + k+pY {q + k + l+ pY{q + lY{q + k+ p) 
6. The same diagrams with the Pauli-Villars loop 



2 



d^k dH dS 



(27r)4 (27r)4 (27r)4 A;2(^1 + (-l)nA;2"/A2n) /2(^1 + (-l)n/2n/A2nj 



X 
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X 



X 



(g2 - Mi) ({q + pY - Mi) ({q + kf - Mi) {{q + k + If - Mi) 

2{q + k + lf[{q + k+p)^ + 



({q + iy - Mi)({q + k+py - Mi)[ 
+ {q + kf-p^)+k\{2q + k + lf + {2q + k + l+pf-p')-mt- 

-2Mi{2q + k + lf -2Mi{2q + k + l+ pf + 2Mi(A{q + k + pf + 

+4(g + A;)' - A;' - p') + 2Mf (5 + l)^,{2q + 2k + 2l- 2p)^ + 

+ 7 r -2q^{q^k^l^pf - 

({q + k+ pf -Mi)({q + k + l+ pY -Mi)[ 

-A{q + pf{q + k + l+pY -8{q + k+ pf{q + pf - Aq^{q + k+pf + 



+2p2 (6Mf - + 2f (2{q + pf - p") + (A; + if {2{q + k+ pf - 



+ 



+2Mi)-6M^ + AMi(2{q + k+pf + 2{q + pf-k^-p^) + 

+ 7 V -QM^ + {2q + k + l+pf X 

({q + l+pf-Mi)[{q + k + l+pf-Mi)[ 

x[-2{q + k + I + pf + {k + If + p"") + AMi{3{q + k + I + pf - p^) - 

-3Mi{k + If - 2Mi{k + qf + 2Mi{l + p)"^ 
16 



+ 



{{q + l+pf-Mi) 
^(q + kf-3Mi)({q + k+pf + 



{{q + kf-Mi){{q + k+pf-Mi) 
+{q + kf + + + pf - e) + 2My + hMf -Mf{q + k) 
2{2{q + k + lf-{k + If) {{q + kf - Mf) AMf 



+ 



+ 



({q + k + pf -Mi)({q + k + l+ pf - Mi) ((g + If - Mi) 
2( - q\q+pf + 6Miq' - 9M/) {{q + k + If - Mj) iQMj 
^ iq' - Mi) ((g +pf-Mi){{q + l+pf- Mi) ^ q' - Mi ^ 
4( - + 3Miq^ - mf) {{q + k + If - Mi) 8Mi{2q + k + If - IGMf 



{q'-Mif{{q + lf-Mi) 
4 



{q2-Mi){{q + kf-Mi) 



[q^-Mi){{q + lf-M. 
-{q + kfq^ + 3Mf -Mi{q + kf - 5M, 
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7. Diagrams containing one insersion of counterterms, presented in Fig. 6: 



/x = -2z - In - + c/i + — In' - + c/iln- -72a'(ln- + c/2 X 

8. Diagrams containing two insersions of counterterms, presented in Fig. 8: 

ia^ ( A n2 
/xx = ^( ln- + ^i) X 



X 



9. Two-loop diagrams containing insersion of counterterms, presented in Fig. 7: 
_ 2v- «/, A \ [ d^l^^ AMI 

I ((A; + qY -Mi){{k + q + pY - Mf) {{q + pY - Mf) (f - Mf) ' 

2(/,- + .y)2 



{{k + qY - Mf) ((g + pY - Mf) (g^ - Mf^ 

4(g + p)2m2 + 2g2M2 - 6Mf 



(79) 



((fc + q)2 _ Mf) ((g + pY - Mf) (f - Mf) J 

B Three-loop contributions to the /^-function. 



In order to find the three-loop /9-function, it is necessary to calculate the integrals, 
presented in Appendix A. 

1. Performing the Wick rotation and using the standard technique, it is easy to 
see, that 

f d^k f 1 ^ 1 

h-ioop-^J ^2^y[j,2^k+pY Y\k^ + Mf){{k + pY + Ml 



1 V- /, 



AMf , / \ I A 



= E ( I. V + ^ 1 + l/-ta.h^ J . _ to _ + 0(1). (80) 

(We assume, that Mj = a^A, where are some constants.) 
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2. Next we analyze graphs containing insersions of counterterms on matter lines: 
Integrals in fx, fxx and fx2 are functions of p/A finite at A — s> oo. Therefore, the 
divergent part of the effective action is defined by their values at p = 0. In this limit 
expressions (77) and (78) in Euclidean space can be written as 




(82) 

where we take into account, that X^Q = 1. In order to calculate fx2 we note, that at 



fa A \ d f d'^k d^q 



X 



I ((A; + g)2 + M2)'(g2 + M^)' ((A; + q^ + Mf) (f + uff ) 
In appendix C we prove, that this integral is equal to 0, and therefore 

fx2 = 0. (84) 

Expressions (81), (82) and (84) agree with the exact result for the sum of diagrams 
with insersions of counterterms, obtained in [30] : 

- In Z Re / d'^x d^O WaC^Wi, + finite terms. (85) 

And indeed, for the considered theory Z is given by (40) and the terms of the considered 
order in a in (85) are 

— - ln- + ^i) + — (ln2-+yiln-) - 72 In - + ^2) - 



/, A 



= {fx + fxx + fx2 + 0{a')) Re J d^x dH WaC^'W,. (86) 
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It means, that the result for the sum of Feynman diagrams agrees with the exact result 
(66). 



3. In Euclidean space foo is given by 



2 f d^k d'^q 



(k+p + qf + q'^ -k'^ -p^ 



X 



X 



+ fc2"/A2"j (k + qf{k + P + qYq^{q + pf 
1 



dH 



2e' 



(27r)4 t^{k + ty 



(27r)4 (t^ + Mf)[{k + ty + Mj 



(87) 



In order to prove that this expression contains only the first degree of In A, it is necessary 
to verify the existance of the limit 



A 



dfoo 



dA 



2 ^ d f d'^k d^q 



^ 2e A— - 

dkJ (27r)4(27r)4^2(^i + ^2n/A2n)'(;t + g)V 



X 



(l + ^-/A-)fl + ^(ln^ + 60)-2e^ 



dH 



47r2 V ^ 
1 



(27r)4 [t^j^Mf)({k + tf + Mf), 
Taking into account, that 



/A 



(l^^{k^, + g/.) _ vr^ 



,2 

( 

dH 



X 



1 



(27r)4 t^{k + ty 



{k + qYq^ P 
(this identity is derived in Appendix C), (88) can be written as 



(89) 



A 



dfoo 



dA 

e 



p=0 



d^k 



87r2 dAJ (27r)4^4(^l+^2n/A2„^ 

dH 1 



-2e" 



r a I 
J (2^ 



(27r)4 + 



(27r)4 (i2 + M2)((fc + i)2 + M2), 



(90) 



It is important to note, that there are some graphs, containing an internal loop or 
insersions of counterterms on the photon line (first 5 diagrams in Fig. 9), contributing 
to the two-loop two-point Green function of the matter superfield. According to [29] 
their contribution in Euclidean space is 
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An d^P 



{2ny 



X 



X 



(l + ,-/A-)(l + £l(ln^ + 60)- 



dH 



Y t\k + tY 



(27r)4 [t2 + Mf)[{k + tY + Mf). 



(91) 



and contains only the first degree of In A. By comparing (90) and (91) we find that, 
the hmit (90) exists and is equal to — I/IGtt^ multiplied by the corresponding two- 
loop contribution to the anomalous dimension. This, in turn, yields the following 
contribution to the /3-function: 



a 



A/3 = A7, 

TT 



(92) 



where A7 is a contribution to the anomalous dimension from a one-loop diagram and 
two-loop diagrams, containing corrections to the photon propagator. 



4. In order to calculate the divergent part of fo we prove, that the limit 



p=0 



(93) 



where fo{A/p) is given by (75), exists. Using the identities 

(k + lf = (q + k + If + q^-{q + kf - (g + 0^ + P + /2; 
{2q + k + lf = {q + k + lf + {q + kf + {q + lf-k''-f + q^ 

in Euchdean space it is possible to present (93) as follows: 



(94) 



X 



A 



d 



r d^l 
J ?2^ 



d^k dH d^q 



2e' 



P=o dA J (27r)4 (27r)4 (27r)4 ^^2(1 + A;2"/A2«) P(l + / A^^) 

( 2 2e ^q + iy 



X 



+ 



q^'iq + lYliQ + kY (Q + kY {q + kfiq + k + lY (q + kfiq + k + lY 
2{k^ + kH'') 2k''{k + lY {k + lY \ 



+ 



+ 



(95) 



{q + kYiq + k + lY {q + kY{q + k + lY {q + k + lYj' 
This expression can be simplified by identities (115) - (117), presented in appendix C 



^dA^^ 



1 



A 



p=Q 
d 



(96) 



d^k dH 



47r2 dA J (27r)4 {2nY (l + /A2») (l + P^i /A2n) [ 2kH^ kH^{k + If 
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This integral is a finite constant (see appendix C). However in order to relate the three- 
loop /9-function with the two-loop anomalous dimension, it is convenient to rewrite (96) 
in the following form: 



2e' 



167r2"dA V J (27r)4 {2tiY kH^{k + p)2(/ + py(l + k^^/K^'^) (l + l^^/k 
d^k dH 4e^ 

+ 



2n 



+ 



(27r)4 (27r)4 ^2^2^^ ^ ^y^^ + / + p)2(l + A;2n/^2n^ ^ ^2ry^2„^ 

_d^_dH 4e\k + / + 2py 

(27r)4 (27r)4 ^^2^^ + p)2/2(/ + p)2(^ + / +p)2(^l + A;2«/A2-) (^1 + /2n/A2n 

because in this case (90) and (97) give 



p=0 

(97) 



p=0 



(98) 



p=0 



where zi and 2^2 are defined by (45) and (46) respectively. 
5- foopv is finite. Indeed, in Euclidean space 



_ sr^ f d'^q 



X 



+ fc2«/A2"y 

{k + p + + q"^ — k'^ — p^ 



X 



k + qf + Mf) [{k+p + qf + M]) {f + M]) ((g + p)2 + M] 

AM] 



X 



((fc + qf + M2) (g2 + Mf) ({q + p)2 + M|) . 



X 



dH 



(99) 



(27r)4 (t2 + ^2j((^ + ^)2 + ^2j^ 

This expression can be written as 

/oop. = /i(p/A) + /2(p/A)ln-, 

A* 

where the functions /i and /2 can be easily found from (99). In appendix C we prove, 
that 



(100) 



r d^q / {k + qf + q^ - k 



(27^ y \ ({k + qf + M2) (^2 + M"^) {{k + g)2 + M^) (f + M2)^ 



0. (101) 
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and therefore /i(0) = /2(0) = 0. Since the functions fi and /2 are evidently holomor- 
phic at = 0, this means that 

hm foo,y = 0. (102) 

A— >oo 

6. The finiteness of fopv can be proven similarly. Indeed, it is evident, that 

f0,y^f{p/^). (103) 

However, at p = identities (118) - (123), presented in appendix C, give /(O) = 0. 
Therefore, fopv is finite and it vanishes in the limit of the regularization removed. 

C Calculation of three-loop integrals, regularized 
by higher derivatives. 



The integral 

can be calculated in four-dimensional spherical coordinates {q,9i,92,(fi)- In these co- 
ordinates we have 



loo = dQ dq -7 

(q^ + 2qkcosa + k'^] 



oc 2 

I dVL I dq---, ^ -^-fdQ^^'^, (105) 

27 7 ^ dq(^q2 + 2qkcosa + k^) '2 J k^ k^' ^ ' 

where a denotes an angle between four-vectors k and q. This angle can be chosen equal 
to 6*1, while 

dn = ddi ddi dip sin^ sin 9^. (106) 
An integral in (83) and (99) at p = 



1X2 



-J (27r)H' 



((A; + g)2 + M2) (g2 + M2) ((fc + g)2 + M2)(g2 + M2) J 

(107) 

can be computed similarly: In the four-dimensional spherical coordinates 
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CXD , 



(27r) 



1 °° 



4M2g3 ] 



0. (10^ 



In order to compute integral (96) we take into account its symmetry with respect 
to the substitution k <-> I and present (96) in the following form: 



p=0 

I . d f d^k dH 



e + p-{k + iy 



(27r)4 (27r)4 (^1 + k^^'/A^^') (l + P^/A^^) 2kH\k + If 



n r d'^k d^l 



g4 /2rY^2n 



7r2 y (27r)4 (27r)4 ^ ^2n/A2n) (1 + ;2n/A2n)^ kH^{k + If 



(109) 



The integral over d^k can be calculated in the four-dimensional spherical coordinates 
if the fourth axis is directed along Z^: 



J j47 1 °° 27r 

d n f d I \ r t t t 



(27r)4 (27r 

e^/2n/A2n 



X 





COS 9i 



(1 + A;2'^ /A2'^) (^1 + /2n /A2n) ^ (A;2 + 2kl COS 9i + /2) ■ 



(110) 



After the substitution x — cos 9i the integral over angles is reduced to the integral over 
contour C, presented in Fig. 10: 



47r / dx 



X 



x^ 



A;2 + 2klx + P 



27r f dx 



xy/1 — ^2 



/ k^ + 2klx + P 
c 



i Resl 



X 



vT 



x^ 



\k^ + 2klx + I- 



-, X = 00 - 



-An'^i Res 



X 



VT 



x^ 



AnH - 



Akl 



+ 



k'^ + 2klx + P' 

ii^ + pf i\e-p\{e + p)\ 



X — 



e + P 
2kl 



8k^P 



8kH^ 



(111) 
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Therefore 



and 



27r 



9 dx ■ 



X 



x^ 



J P + 2klx + P 
c 



k>l; 



I > k 



(112) 



ne 



dl 



xi j dk 



X 



+ / dA; 



A; 



k^(l + A;2" /A2») i ^3 (^1 + ^2n /A2n) 



Making the substitutions a; = /^/A^; y = A^/A;^ in the first integral and x 
y = /c^/A^ in the second one, we finally obtain 



(113) 

AV/2; 



A ^ , 

^dK^- 



1/x 



ne 



dx 



p=o 327r2(27r)4y (1 + a;")^ J 1 + y-" 



ne 



oc 



ne 



327r2(27r) 



dx 



X 



n-l 



(i + x'^y 



327r2(27r) 



327r2(27r) 



dx 



X 



l/x 

- dy 



I 



(1 + x")2 y 1 + 



327r4' 



(114) 



Below we also present the identities, which were used for taking integrals, which 
contain higher derivatives. 

1. Identities required for calculation of diagrams with 0- and 0- lines: 



Ii = 



d^k dH d^q 



{k + qf + q^ -k"^ 



(27r)4 (27r)4 {2tiY k^[l + A;2"/A2n) /2(l + /2n/A2n) qi^q + kY{q + If 
1 r d^k dH 1 



167r2 j (27r)4 {2tiY (^i + A;2'^/A2'^) (l + /2n/A2n) kH''' 



(115) 



d^k dH d'^q 



(27r)4 (27r)4 (27r)4 ^^2(1 + A;2n/A2") P(l + / A?-^)q\q + Z)^ 
12A;2 2(A;4 + A;2Z2) 2k\k + lf 



X 



(g + A;)2(g + A; + 0' (g + A;)4(g + A; + Z)' (g + A;)2(5 + A; + 0' 



27 



1 r d^k dH 1 

4^ J (27r)4 (27r)4 (l + k^n//^2n^ (l + /2n/A2nj ^4^2 + ;)2' 



"i (27r)4 (27r)4 (27r)4 A;2(i + A;2"/A2") ^2^1 + Z2n/A2„^ ^4(^ + ^ + ;)4(^ + 

1 r d^jH 1 [ 2 ]_1 

167r2 7 (27r)4 {2nY (l + /A2n) (1 + /A2n^ | kH'^{k + Z)2 ~ 'PI? J 



(117) 

2. Identities, required for calculation of diagrams with internal Pauli-Villars lines: 



J _ r d^k dH dS 1 ^ 

' " y (27r)4 (27r)4 (27r)4 ^^2(1 + k^n^/^2n^ 12(1 + i2n/A2nj + ;)2 + 7^4-2 j ^ 

X — ——T-. — — :v + 



{q^ + M2)2 ((5 + kr + M2) (^2 + ^2) ((^ + ^)2 + M^)' 
k'^ + 2M2 2M2 



(g2 + M2)2 ((g + A;)2 + M2) ' (?2 + M2)3 ((5 + A;)2 + M^) 

d'^k dH dS 1 



0; (118) 



J2 = 



(27r)4 (27r)4 (27r)4^2(^i + k^^/A^^) P(l + P'^ / K^''){q^ + M^f 
1 f k\e + P) 



((g + A;)2 + M2) ((g + Z)2 + [ ((g + A;)2 + M^) ({q + k + If + M2) 

2M2(A';2 + /2) 4M2A;2 



((g + A;)2 + M2) ((g + A; + 0^ + M^) {q^ + M^) ((g + A; + If + M2) 
6A;2 k\k + lY + 2MV \ ^ 



+ 



((g + A; + 0^ + M2) ((g + A; + 0^ + M'^)' 



(119) 



/■ (i^fc rf^g 1 

y — / ^ 

^ ~ J {2nY (27r) i (27r)4 + A;2"/A2") /2(^i + ;2„/A2n)(g2 + M2)2 

1 f (A; + /)2 + 2M2 4M2 



((g + A;)2 + M2) ((g + A; + /)2 + M2) I [[q k + If + M^) {q^ + M^) 
{q + k + iy + M2 A;2 + 2M2 



((g + /)2 + M2) ((g + A;)2 + M^) 



0; (120) 
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Ja 



d^k dH dS 



(27r)4 (27r)4 (27r)4 ^2(^1 + k^n//^2n'^ 12(1 + /2n/A2n)(^2 + JVf2)^ 



1 r^,„ M\k + lf + 2M^ 

X7 TTT -<^4M^ - 7 ^ v + 

((g + ky + M2) {{q + k + ly + M2) I ((? + A; + Z)^ + M^j 

((g + A;)2 + M2) ((^ + A;)2 + M^) (g^ + M2) j 



, d'^q 1 
J. = / — 7 ^ 7 ^ X 



/ 



{2txY (27r)4 (27r)4 k^l^l + fc2"/A2'^) /2(^1 + /2n/A2n) ((g + /)2 + M^) 

2M2 M2 



(g2 + M2) ((g + A;)2 + M2) (g2 + M2) ((g + A;)2 + M2)' 
2M^ + M2A;2 SM^ 1 



(g2 + M2) ((5 + A;)2 + M2) (g2 + M2) ((g + A;)2 + M2) 



0; (122) 



f d'^k dH d^q 1 
^ ~ y (27r)4 (27r)4 (27r)4 ^^2(1 + A;2-/A2-) p(l + /2n/A2n^ (^2 + 7^2)2 

1 f 2M^ 

((g + fc)2 + M2) ((g + /)2 + M2) 1 (g2 + M^) {{q + k + If + M") 

M\{q + If + M2) - M2 2M2 



((g + A; + 0' + M2)' ((g + + 0' + ^2)' [{q + k + lf + M^) 



0. (123) 



((g + kf + M2) ((g + A; + If + M2) , 
As an example wc prove identity (115). For this purpose we consider 

^ ^ rJ^J^^ 2(g^ + qM 

' J {2Trf (27r)4 {2Trf k"" (l + A;2"/A2«) P (l + /2n/A2n) g4(^ + kfiq + lf 

and write the integral over d'^q in four- dimensional spherical coordinates: 



d^k dH 1 

X 



(27r)4 (27r)4 ^^2(^1 + A:2«/A2«) /2(^i + /2„/a2 

X A [d^ hi ^{<l + <lk CO. a) 

[2nfJ J (g2 + 2gA;cosQ; + A;2) (g2 + 2gi cos /? + ^2) 

d^k dH 1 

X 



(27r)4 (27r)4 ^^2(^1 + A;2"/A2") ^2(^1 + Z2n/A2n^ 
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Iff 1 d 1 

/ dn / dq-, ^ — 7 r-(125) 

(27r)4y J {q^ + 2qlcosP + P)dq[q^ + 2qkcosa + k^) 

Performing integrating by parts in the last integral, we obtain 
d^k dH 1 



(27r)4 (27r)4 ^^2(^1 + k^^/A"^^) P(l + /2n/A2n^ 



X 



X— ^ fdQ ^ 



(27r)4 y I (g2 + 2qk cos a + A;^) (g2 + 2g/ cos /3 + p) 



(g2 + 25A; cos a + k'^) dq (g2 ^ 2ql cos /3 + Z^) 

r d^k dH 1 

J (27r)4 (27r)4 ^^2(1 + k^'^/A'^^) P(l + / K^^) ^ 



X 



1 f 1 °° 



2(g + A;cos/3) 



(27r)4y [A;2;2 j (g2 + 2gA;cosa + A;2) (g2 + 2g/cos/3 + /2)^ 

d^k dH 1 
(27r)4 (27r)4 k^l^i + k^nj^in^ 12(^1 + /2n/A: 



/ 



l87r2A;2/2 7 (27r)4 g4(g + A;)2(g + /)4 f ^ > 



The last integral in this expression is evidently equal to /i. Therefore, 

' 167r2 J (27r)4 (27r)4 (1 + k^^n j p,2n^^ (1 + /2n/A2n^ ^4^4 " ^ ' 

The other identities can be derived in the similar way. 
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Figure 1: One- loop graphs for calculating AFy. 




Figure 2: Two- loop graphs for calculating AFy. 




Figure 3: Three- loop graphs with two loops of the matter fields. In these diagrams 
external lines are attached to both loops. 
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Figure 6: Diagrams, containing an insersion of two- loop counterterms on the matter 
lines. 
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